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1 Abstract 
Os 

0$ ■ We study the topological defects in the classical Heisenberg antiferromagnet 

S 

in two dimensions on a triangular lattice (HAFT). While the topological anal- 



ysis of the order parameter space indicates that the defects are of Z2 type, 
consideration of the energy leads us to a description of the low-energy sta- 



^ ■ tionary points of the action in terms of ± vortices, as in the planar XY model. 

Starting with the continuum description of the HAFT, we show analytically 
that its partition function can be reduced to that of a 2-dimensional Coulomb 
gas with logarithmic interaction. Thus, at low temperatures, the correlation 
length is determined by the spinwaves, while at higher temperatures we ex- 
pect a crossover to a Kosterlitz-Thouless type behaviour. The results of recent 
Monte Carlo calculations of the correlation length are consistent with such a 
crossover. 

PACS numbers: 75.10,75.10. J, 75.50.E 
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The effect of frustration on the properties of two-dimensional antiferromagnets has been 
the subject of numerous investigations [1-7]. The Heisenberg antiferromagnet on a triangular 
lattice (HAFT) is a prototypical model for such systems. In early studies the focus has 
been on the ground state properties of the quantum mechanical model. These properties are 
qualitatively similar to those of the same model on a square lattice (HAFSQ) ||[| where 
there is no frustration. What will be of interest here and possibly also for the interpretation 
of experimental results is the behaviour at finite temperatures. As is frequently observed, the 
HAFT possesses vortex excitations, which are important at finite temperatures and which 
are absent in the case of the HAFSQ. The reason is the difference in the order parameters. 
On the triangular lattice, in contrast to the square lattice case, the order of the spins at 
zero temperature is not collinear; instead, the frustration leads to a state in which, on every 
elementary triangular plaquette, the spins form a regular planar star. As a consequence, 
the order parameter space is isomorphic to the full 3-dimensional rotation group SO (3) 
(HAFSQ: SO(3)/SO(2) = S2). The 5*0(3) (in contrast to the S 2 ) is not simply connected, 
and a homotopy analysis predicts the existence of topological stable point defects. Because 
the fundamental homotopy group is n 1 (S'0(3)) = Z 2 [|TU[ these defects are called Z 2 -vortices. 
The earliest hint as to the importance of the vortices came from the Monte Carlo simulations 
of the classical HAFT by Kawamura and Miyashita ||. They found evidence for a defect- 
driven transition which is strongly reminiscent of the Kosterlitz-Thouless transition [|ll|] in 



the planar XY model. However, there, in the 5*0(2) case, the vortices are characterized by 
integer numbers (n 1 (S'0(2)) = Z) ||10|| . From their data, Kawamura and Miyashita cannot 
conclude whether the transition is of Kosterlitz-Thouless type or if it is caused by a different 
mechanism. 

The purpose of this letter is to treat analytically the vortices in the HAFT. It will be 
shown that, after an integration of the spinwaves in harmonic approximation, the partition 
function reduces to that of a 2-dimensional Coulomb gas with logarithmic interaction. Thus, 
we arrive at a similar description as in the case of the XY model. In the derivation of this 
result we have to face two problems, (i) How can one understand the picture of a gas of ± 



charged particles, when one knows from homotopy considerations |10[ that there is only one 
type of vortices ? (ii) In constrast to the case of XY model, the spinwaves in the HAFT 
are coupled to the vortices, already in the harmonic approximation. While the vortices 
interact logarithmically in the absence of this coupling, it is not at all trivial to show that 
this logarithmic form survives despite of the coupling. Furthermore, anharmonic spinwave 
interactions are known to yield a finite correlation length £$w for the HAFT for arbitrarily 
low temperatures whereas the correlation length of the XY model is infinite below the 
Kosterlitz-Thouless temperature. As we shall discuss below, this restricts the validity of the 
logarithmic vortex interaction to distances r <C £,sw- 

Analogous to the investigation of the planar XY model, our starting point is the action 
of the continuum version of the classical HAFT, i. e. the nlcx model [p,[12 



S{n k } = -lJd^{(Vn 1 ) 2 + (Vn 2 ) 2 } , t = ^% . (1) 

The orthonormal vectors ni(x) and n^ix) span the plane of the original planar star of 
spins on each elementary triangular plaquette. Since the 577(2), as the universal covering 
group of the 50(3), plays an essential role in homotopy analysis |T0|] , we represent the field 
configuration by SU(2) fields g(x): 

n k = l -tr{aga k g- 1 } , g(x)eSU(2) . (2) 

Here a k (k = 1, 2, 3) are the Pauli matrices and % = nxxn 2 . In general, the SU(2) fields can 
be written as g(x) = exp || u(x)-a \l/(x)| where u(x) is a unit vector specifying a local axis 
of rotation and \l/(x) is a local rotation angle. From (^), one can see how a configuration 
nfc(x) is generated by local rotations g(x). The ambiguity in the sign of the SU{2) fields in 
(0) reflects the existence of two classes of states of different topology. 

Next, we determine the stationary points of (^). The factorisation of the fields 
g(x) = g s (x) exp{i e(x) ■ a} into a saddle point and a spinwave part leads to Euler-Lagrange 
equations. The trivial solution g s (x) = const reproduces the spinwave minimum. Searching 
for nontrivial solutions g s {x), we assume u(x) = u s = const: g s (x) = exp|| u s -a ^ s (x)\. 
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Minimal action solutions are found to be those which satisfy u s -e z = 0. This ansatz reduces 
the Euler-Lagrange equations to the Laplace equation for the field *& s (x). Thus, ^ s (x) 
must be identical to the vortex configurations of the SO(2) model. The one-vortex solution 
centered at the origin is: ^i(x) = arctan a ; g\{x) = exp{| u s - a arctan a }. We have proven 
local stability against small fluctuations e(x) explicitly. The energy of one vortex diverges 
logarithmically with the system size. 

From the topological point of view, a Z2 vortex can be characterized by the behaviour 
of the field g(x) on an arbitrary path encircling the singularity: as one moves around such a 
closed path, g(x) picks up a phase factor e l7T , still the fields are continuous corresponding 
to the above mentioned ambiguity in the sign of the £77(2) fields in (0). A general one- 
vortex configuration is obtained by a continuous deformation of the stationary configuration 
g\{x). However, in contrast to the £0(2) case, the stationary point g\ in the homotopy class 
with one vortex is degenerate with respect to the direction of u s . Solutions with different 
orientations u s can be continuously transformed into each other. In the £0(2) case, vortices 
and antivortices are stationary points of the action in different homotopy classes. 

Since the Laplace equation for ^f s (x) is linear, stationary N-vortex configurations are 
obtained simply by linear combination: 

^n{x) = ^2 s j arctan^ — — , Sj = ±1 (3) 
3=1 x ~ x j 

with A singularities at x±, ...,xn and (?iv(x) = exp || u s -a tyjy(x)\. Finiteness of the action 
demands YljLi s j — 0, i. e. has to be even. The value of the action at the stationary point 
is: 

4 0) = ft E W lu^^ , (|x,-x>| »a) . (4) 

3,3' 

The vortices interact logarithmically (a is the lattice spacing). As in the one-vortex case, 
the A-vortex stationary point is degenerate with respect to the direction u s . For any 
fixed u s , such a configuration consists of an equal number of vortices (rotation axes u s ) 
and antivortices (rotation axes —u s ). A more general configuration with an even number 
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of vortices, g^, is obtained by forming a product of N single vortex fields gi(x — Xj) = 
exp{| u s j- a arctanl^f 2 "} {j — !)••■> -^0 with independent directions u s j of the individual 
vortices. However, in general, this is not a stationary point of (|I|). We have checked for 
small iV (< 4) that for arbitrary directions u s j the energy of such configurations is always 
higher than (^). Thus in addition to homotopy theory, one needs to consider the energy, 
which then singles out the configurations of the form of with equal number of vortices 
and antivortices of a fixed axis u s . This is the solution to the first problem posed in the 
introduction. 

We now approximate the partition function of our model by summing over the stationary 
points and including small fluctuations up to second order: 



N=0 If! 



Z ~ fN \2 j d 2 x 1 - ■ ■ j d 2 x N Z N (x ir ■ ■ ,x N ) (5) 



Z N (x h ■ ■ ■ ,x N ) := expSP jv[e\ exp|«S^ ) [e(5)]| , (6) 
where z is the fugacity of a pair of vortices. The fluctuation part of the action reads 

S% ] = ~ Jd 2 xe(x) T M N [x\ e{x) . (7) 
In the basis (u s , e±, e z ), M N [x] is the 3x3 differential operator 



M N [x\ :-- 



' -A ^ 

-A + (W w ) 2 -2(Wjv) V 



(8) 



i 2 (Wat] V -2A j 



The first component of the fluctuations generates small local rotations around the axis u s . It 
corresponds to the spinwaves of the 5*0(2) model and decouples from the vortices as there. 
The other two components couple to the vortices. Without this coupling, the vortex part 
of the partition function would be identical to the vortex part which one gets for the planar 
XY model. A priori, it is not clear, how the fluctuations influence the thermodynamics of 
the vortices. In order to study this second problem posed in the introduction, we need to 



derive an effective description of the thermodynamics of the vortices. To this end, we have to 
integrate the fluctuations in eqn. (|5]) in the presence of the vortices, i. e. we have to calculate 
the functional determinant of the semi positive elliptic operator Mn[%\- This is a formidable 
problem, and to our knowledge, such a problem has been solved only in a few cases. For 
example, Fateev et. al. calculated such functional determinants in their investigation of the 



thermodynamics of the instantons in the SO (3) nkr model [13 . 

Now, we first concentrate on the case of two vortices, N = 2. The contribution of 
the fluctuations to the effective interaction potential of these two vortices be \nW : = 
ln{detM2/detMo}. In order to determine that part of the functional determinants which 
diverges with the number of degrees of freedom, we used the proper-time regularization 



and calculated the first nontrivial Seeley coefficients via a Heat-Kernel expansion |TJ[]. The 



leading divergence cancels in numerator and denominator of W. The next to leading (loga- 
rithmic) divergence vanishes exactly, i. e. W is finite. 

We further evaluated W, using two different regularizations. In the first, the analytical 
calculation, we observe that in bipolar coordinates, the differential operators in M 2 separate. 
The bipolar coordinates are ip and p []15|| , where <p(x) = ^f%(x) and p(x) is the orthogonal 
coordinate. The x-y plane maps onto a cylinder of length 21n^ where r = \x± — X2\- 
Furthermore, the coefficients of the differential operators in M2 become constant. Gaussian 
integration of the eigenmodes of nonzero energy and exact integration of the Goldstone 
modes over the group space of SO(3)xSO(2) yields an expression for \iaW containing sums 
over logarithms of the eigenvalues. Extracting the asymptotic for r>awe finally obtain 

lnWgr« W ln£ + ^lnln£ + 0(l) • (9) 

The expression for w evaluates to 0.06. This result shows that, within the harmonic approx- 
imation (0), the coupling of the fluctuations to the vortices does not change the asymptotic 
logarithmic behaviour of the interaction of the vortices. 

We confirmed this result in a second regularization by calculating In W numerically on 
finite lattices in real space. We computed In W for lattices up to sizes of 64 x 64. For the 
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values of r (r < 64) available, the results converge towards the corresponding result of the 
analytic diagonalisation. But, owing to the small value of w, one cannot reliably identify 
the asymptotic In - behaviour in the numerical data. To obtain further confirmation of the 
agreement between the analytic and the numerical calculation, we generalized the action ([!]) 
by adding a term ai^Jn?) 2 which does not change the symmetry (see ||). The calculations 
were then also performed for several values of a with < a < 1, where a = 1 is a system of 
the higher symmetry 5*0(3) x 5*0(3). Increase of a from to 1 increases the prefactor w in 
(§) smoothly to values of the order of unity. This makes it possible to identify for a > the 
In - behaviour in the numerical data even for small values of r. Especially the dependence of 
w on a shows a very good agreement between the analytic and the numerical calculations. 
However, the magnitude of the prefactor of the leading term in @ is less important. The 
essential point is that, within the harmonic approximation, the fluctuations do not change 
the logarithmic form of the interaction potential of two vortices at large distances. 

Finally, we turn to the case of N > 2 vortices. As in the 50(2) model, the two-vortex 
excitations are the elementary local excitations. At large distances, the interaction between 
such excitations decays as the inverse square of their distance and should not be changed 
by the fluctuations. With this assumption, we get for the general N-vortex case 



ln ^(f 1 ,---,f JV ) = 2 ! r J^-xA 



([10|) together with @ leads to an expression for the partition function which describes a 



2-dimensional Coulomb gas ||11|| . The connection between the dimensionless temperature t 
and the effective temperature of the Coulomb gas t' is: h = ■h (1 — f-t). 



The representation (^) with ( |TOD for the thermodynamics of the vortices constitutes the 
result of our work. In the following summary, we want to describe the implications for 
the classical HAFT. For small t, the Coulomb gas (|T0| ) is in its low temperature phase. 
The Kosterlitz-Thouless transition temperature is t' KT = 2.7. Thus, one would expect a 
transition in the HAFT at around txr — 0.68 (or Tkt — 0.58J5 2 ). However, in the steps 
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leading to (|10|), anharmonic spinwave interaction have been neglected, i. e. we considered 
a system of vortices in an ideally correlated background. As we have pointed out above, 
these anharmonicities lead to a finite correlation length £$w even in the absence of vortices. 
Hence, the validity of our results is restricted to distances smaller than £,sw- This, of course, 
implies that at T KT the system does not become critical. Nevertheless, assuming that ^sw 
is sufficiently large, we expect a spontaneous creation and dissociation of pairs of vortices 
to occur as the temperature increases towards Tkt, resulting in a sudden breakdown of 
the remaining correlations. The temperature dependence of the correlation length should 
then cross over from the low temperature behaviour £sw ~ Viexp(l + \)ir/t to the 



Kosterlitz-Thouless behaviour £kt ~ exp b/ \/t — txr- Regarding the difference in the value 
of our estimate of Tkt and that of MC simulations of Kawamura and Miyashita ( ||: 
Tkm — 0.31 iO.OUS 12 ) we refer to the similar situation in the case of the planar XY model. 
Tobochnik and Chester found in their MC simulations of the planar XY model |16|] a value 
for the transition temperature Tkt which is about 30% reduced from the value which was 
found from Kosterlitz and Thouless. 

It is interesting to note that at tKT the spinwave correlation length £sw would be very 
large (= 10 5 lattice spacings) if one assumes the unknown prefactor to be of the order of 
unity as for the classical HAFSQ [[!]]]. In recent MC calculations of the classical HAFT, 
Southern and Young Jl8j] found a sharp drop of the correlation length at a temperature 
which agrees with the estimate Tkm of Kawamura and Miyashita. For T > Tkm, the data 



in 



TEH can De fitted well using £kt- The interpretation of this behaviour in ref. [|lj| is the 
same as presented here: there is a crossover in the correlation length of the HAFT due to 
the unbinding of the vortices. 
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